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In this paper we demonstrate how, using a natural generalization of BCS theory, superconducting phase coher- 
ence manifests itself in phase insensitive measurements — when there is a smooth evolution of the excitation 
gap A from above to below Tc- In this context, we address the underdoped cuprates. Our premise is that just 
as Fermi liquid theory is failing above Tc, BCS theory is failing below. The order parameter Asc is different 
from the excitation gap A. Equivalently there is a (pseudo)gap in the excitation spectrum above Tc which is also 
present in the underlying normal state of the superconducting phase. A central emphasis of our paper is that the 
latter gap is most directly inferred from specific heat and vortex core experiments. At the same time there are 
indications that fermionic quasiparticles exist below Tc so that many features of BCS theory are clearly present. 
A natural reconciliation of these observations is to modify BCS theory slightly without abandoning it altogether. 
Here we review such a modification based on a BCS-like ground state wavefunction. A central parameter of our 
extended BCS theory is A^ — A^^. which is a measure of the number of bosonic pair excitations which have a 
non-zero net momentum. These bosons are present in addition to the usual fermionic quasiparticles. Applying 
this theory we find that the Bose condensation of Cooper pairs, which is reflected in Asc, leads to sharp peaks 
in the spectral function once T <Tc. These are manifested in ARPES spectra as well as in specific heat jumps, 
which become more like the behavior in a A transition as the pseudogap develops. We end with a discussion of 
tunneling experiments and condensation energy issues. Comparison between theoretical and experimental plots 
of Cv, and of tunneling and vortex core spectroscopy measurements is good. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.25.Bt, 74.25.Fy 
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I. INTRODUCTION 

In the underdoped regime of high temperature supercon- 
ductors it is now clear that Fermi liquid theory is failing 
and the "smoking gun" for this failure is a (pseudo)gap 
in feaiiionic excitation spectrum above Tc- Many would 
argueEli3 that this failure is evidence for spin-charge sepa- 
ration. However, the fact that the excitation gap evolves 
smoothly into its counterpart in the superconducting phase 
may also be interpreted as evidence for "preformed pairs". 
In this way superconducting pairing correlations are respon- 
sible for the breakdown of the Fermi liquid state. This picture 
appears rather natural in view of the notably short coherence 
length ^ in these materials, which leads to a breakdown of the 
strict mean field theory of BCS. Within this short ^ scheme, 
one considers that pairs form at temperature T* and Bose con- 
dense at lower temperature Tc. These werfiot true "preformed" 
or bound pairs but rather long liveduij pair states. ^Many 
have argufidJor this viewpoint from experimentalistsHTlO to 
theorists yy p j-. 

Our contributiontlTlZl to this body of work has been to 
show how to microscopically implement this preformed pair 
approach at all temperatures T < Tc, by deriving an extep^ 
sion of BCS theory, based, on the ground state of Leggett.Efl 
We have also addressedEm the behavior above Tc- In this ex- 
tended BCS approach, the fermionic excitation gap A (which 
evolves smoothly from above to below TcE3E3) and the mean 
field order parameter Asc (which is non-vanishing below Tc) 



are not necessarily the same at any non-zero temperature; this 
is a reflection of the distinction between T* and Tc- Since 
A ^ Asc, we say that there are pseudogap effects below Tc- 
The normal state underlying the superconducting phase is not 
a Fermi liquid. The excitations of the system can be viewed as 
a "soup" of fermions and pairs of fermions (bosons). The lat- 
ter are very long-lived at and below Tc in the long wave-length 
limit; their number is associated with the difference A^ — A^^. 

This background sets the stage for the important questions 
which we address in this paper. What are the signatures of 
Tc, in thermodynamical quantities such as the specific heat, 
Cv, given the smooth evolution of the excitation gap? How do 
we understand the abrupt appearance of long lived, fermionic 
"quasiparticles" below Tc, and their implications for the elec- 
tronic spectral function A(k, w)? If the superconducting state 
is not Fermi-liquid based, then how does one extrapolate the 
"normal state" below Tc in order to deduce such thermody- 
namical properties as the condensation energy? 

One of the central observations underlying this paper is 
the fact that there are two distinct experiments which provide 
seemingly similar information about the extrapolated normal 
state, i.e., that this T < Tc state contains an excitation gap. 
These are scanning tunneling microscopy (STM) dal^ in a 
vortex coreffiil as well as specific heat measurements.E3 Ren- 
ner and co-workerg3 have argued that their vortex experiments 
"show either the presence of importantpSuperconducting fluc- 
tuations or preformed pairs". LoramcacJ and co-workers 
have analyzed their specific heat data to show that a ther- 



modynamically consistent picture of Cy reflects a gap in the 
(T < Tc) "normal" state spectrum, not directly related to the 
condensate. This would, they argue, include the possibility 
of preformed pairs that retained their structure below T^, as 
in He^, and whose binding energy does not contribute to the 
condensation energy. Alternative explanations for the vortex 
core experiments have been advanced by Fjianz and Millis,E3 
and ruore recently by Franz and TesanovicEj and by Lee and 
Wen.B 

This "preformed pair" picture, which is essentially a mean 
field-based approach, should be contrasted, with the phase 
fluctuation picture of Emery and_Kivelsonc3 which has been 
implemented by Franz and MillisEd to address photoemission 
and tunneling spectroscopies. In the present case the normal 
state excitations represent a "soup" of fermions and bosons, 
whereas in the approach of Ref. 23 the system is thought to 
consist of a soup of fluctuating vortices. It is widely believed 
that, at least below Tc, fermionic quasiparticles are present 
so that the phase fluctuation picture wilLeventually need to 
accommodate their contributions. LoramEahas, moreover, ar- 
gued that phase fluctuations are not consistent with the behav- 
ior of C„, which he observes. 

The results which we obtain in this paper show that upon 
entering the superconducting phase, the onset of the coher- 
ent condensate, characterized by lS.sc, leads to a sharpening 
of the peaks in the electronic spectral function, which will 
be directly reflected in ARPES studies |Where its effects are 
quite dramatic, as well as in tunneling.EB ARPES measure- 
ments support such a peak sharpening, and it has been recently 
claimed,E3 as is consistent with the theme of this paper, that 
the observed sharpening at Tc (rather than at T*) is difficult to 
understand within strict BCS theory. 

Indeed, BCS theory can and should be generalized, and in 
its more general form, this peak sharpening in conjunction 
with the temperature dependence of the excitation gap, is also 
responsible for a specific heat jump. The latter, is thus, quite 
generally, associated with the onset of off-diagonal long range 
order When this general picture is applied to the cuprates 
we find that in the overdoped regime, this jump (in Cy) can 
be quantified in terms of the temperature dependence of the 
excitation gap A (as in the traditional BCS case, see Eq. |l8| 
below). This is in contrast to the underdoped regime, where 
the order parameter and excitation gap are distinct and where 
the excitation gap is smooth across Tc- Here the jump (asso- 
ciated only with lS.sc) becomes smaller towards underdoping, 
where the pseudogap is more prominent. Moreover, the shape 
of the Ct, versus T curve is more like the A transition of Bose- 
Einstein condensation (BECXpAll of these features seem to be 
consistent with experiment.c3El 

An important second theme of this paper is an analysis 
of the extrapolated normal state below Tc. We argue here 
that the superconductivity is non-Fermi liquid based and that 
this has important implications for condensation energy esti- 
mates. Moreover, the non-Fermi liquid characteristics of the 
extrapolated normal state (which underlies the superconduct- 
ing phase) should help provide constraints on the long stand- 
ing controversy! of Fermi liquid break-down iBj-the normal 
state. Indeed, spin-charge separation scenarioscZHHa might be 



distinguishable from alternatives such as the present one, or 
that of Ref. ^ by studies below Tc- This provides a major 
impetus for the present work. 



II. THEORETICAL FRAMEWORK 



A. Overview 



Our work begins with the ground state wavefunction of 
BCS as generalized by Leggettta 



*o = nk(uk + t'k4TC-ki)|0) 



(1) 



which describes the continuous evolution between a BCS sys- 
tem, having weak coupling g and large ,f , towards a BEC sys- 
tem with large g and small ^. Here Mk, f k, which are defined 
as in BCS theory, are self consistently determined in conjunc- 
tion with the number constraint. The central approximation 
of this paper is the choice of this ground state .wavefunction. 
The essence of our previous contributionsli3li3'EZl has been a 
characterization of the excitations of ^Pq and their experimen- 
tal signatures (for all T < Tc). New thermodynamical effects 
stemming from bosonic degrees of freedom must necessarily 
enter, as one crosses out of the BCS regime, towards Bose- 
Einstein condensation. 

As in BCS theory, we presume that there exists some at- 
tractive interaction between fermions of unspecified origin 
which is written as Vk.k' — g(pk^k', where g < 0; here, 
ipk — 1 and {coskx — cos ky) for s- and d-wave pairing, 
respectively. The fermions are assumed to have dispersion, 
Ek == 2t||(2 — COS /ca;— COS fcj,) + 2i^(l—cosfcj^) — /i, measured 
with respect to the fermionic chemical potential /i. Here i|| 
and tj_ are the in-plane and out-of -plane hopping integrals, re- 
spectively. In a quasi-two dimensional (2D) system, t± <C i||. 
For brevity, we use a four momentum notation K = (k, iuj), 
J2k = ^X^kw ^^^-^ ^"'^ suppress ipk until the final equa- 
tions. 

We now make a number of important observations about 
BCS theory. BCS theory involves a special form for the 
pair susceptibiUty xiQ) — ^K^i^)^oiQ ~ ^)' where 



the Green's function G satisfies G 



= Gn 



S, with 



I](iir) = —b?gJjQ[—K). In this notation, the gap equation 
is 



1 



5X(0) =0, T < Tc. (2) 

As was first observed by Kadanoff and Martin£il this BCS 
gap equation can be rederived by truncating the equations of 
motion so that only the one (G) and two particle (T) propa- 
gators appeared. Here G depends on E which in tiim depends 
on T. In general T has two additive contributions,Ej from the 
condensate (sc) and the rpn-condensed (pg) pairs. Similarly 
the associated self energyEil 



S(i^)=^r(Q)Go(Q-if) 



(3) 



can be decomposed into T,pg{K)+T,sc{K). The two contribu- 
tions in S come respectively from the condensate, TsdQ) = 
-AlJ{Q)/T, and from the Q ^ pairs, with Tpg{Q) = 

5/(l + 3X(Q)). 

More generally, at larger g, the above equations hold but 
we now include feedback into Eq. (g) from the finite mo- 
mentum pairs, via T.pg{K) = ^Q^pg{Q)G^{Q - K) k. 

Ga{-K)Y.Q%9iQ) = -^pgGoi-K), which defines a 
pseudogap parameter, Apg. This last approximation is valid 
only because (through Eq. (||)), Tpg diverges as Q -^ 0. In 
this way, I]pg{K) has a BCS-like form, as does the total self 
energy E(i^) = -A^Goi-K), where 



A2 = A2 



pg' 



For the physically relevant regime of mO|d^paj;e g 
found, after detailed numerical calculations 
be approximated as 



sBS'thatT;, 



Tpg^{q, n) = ao{Q.-flq+ Upair + i^q) 



(4) 

we have 

pgmay 

(5) 



where the pair dispersion i7q = q^ jIMpair and the effective 
pair chemical potential [ipair = for T < Tc- The effective 
pair mass Mpair rapd the coefficient oq are determined via a 
Taylor expansiorO of X,~^. Moreover, Fq ^ 0, as q ^ 0. 
As a consequence we have 



A^ 
pg 






(6) 



q#0 



We now rewrite Eq. 
straint, as 



b, along with the fermion number con- 



1 



52^ TTB "^k 



E 



k 

Sk 



2^k 
2ek 



Eu 



fiEw) 



(7) 
(8) 



Here f{x) and b{x) are the Fermi and Bose functions and 



9^ is the quasiparticle dispersion. 
) represent the central equations of our the- 



Equations (|3)-(^ 
ory below Tc- They are consistent with BCS theory at small 
g, and with the ground state ^q at all g; in both cases the right 
hand side of Eq. (ph is zero. The simplest physical interpreta- 
tion of the present decoupling scheme is that it goes beyond 
the standard BCS mean field treatment of the single particles 
(which also acquire a self energy from the finite q pairs), but 
it treats the pairs at a self consistent mean field level. 



III. SPECTRAL FUNCTIONS, DENSITIES OF STATES AND 
SPECIFIC HEAT 

Experimentally, it has been established from specific heat 
measurements in the cuprates that there is a step discmvtiriuity 
or a maximum at T^, depending on the doping level.EjO It is 
clear that one cannot explain these experiments using the stan- 
dard picture of BCS theory, in which the specific heat jump at 



Tc results from the opening of the excitation gap. We now 
address these experiments. In Section II, we used an approx- 
imate form for the pseudogap self-energy Y.pg (see derivation 
of Eq. (^), in order to simplify the calculations. Under this 
approximation, Upg has a BCS-like character, so that the spec- 
tral function is given by two (5-functions at ±i?k- These ap- 
proximations were iustified in the context of the applications 
considered, thus far.EjilZl 

However, in order to study quantities which rely on details 
of the density of states, we will, in the remainder of this paper, 
relax this simplifying approximation and allow for lifetime ef- 
fects in Spg. This more realistic form for Spg incorporates a 
finite broadening 7 due to the incoherent nature of the finite 
center-of-mass momentum pair excitations. In this way we 
distinguish this contribution from thatpof the condensate. To 
make numerical calculations tractable,QiJ we do not solve for 
the broadening 7 and excitation gap A self-consistently. This 
would involve an iterative solution of the complex set of three 
coupled equations for G and T; rather|-|We take A below Tc 
from Eqs. (ra)- (M), and use our estimatestacl of T* to determine 
where pseudogap effects are essentially negligible. More im- 
portantly, we treat the broadening as a phenomenological pa- 
rameter, with one adjustable coefficient, chosen to optimize 
fits to tunneling experiments. Our results, throughout this pa- 
per, are not particularly sensitive to the detailed fornjupf 7 
(which will depend on doping concentration x and T)J23 but 
it is essential that 7 be non-zero and appreciable when com- 
pared with T. In this same spirit, we take Tc, and the chemical 
potential /i from our leading order calculations (with 7=0). 

We turn now to the spectral function A(k, Lj).p|U follows 



from our micrO|Sc.opic scheme that slightly aboveaEl Tc, and 
for all T < Tcf3 the self energy associated with the q 7^ 
pairs, and that from the condensate or q = pairs is given 
respectively by 



I]pg(k,w) 



k,pg 



ek + »7 



-«So(k,tj) 



and 



Ssc(k,cj) 



Ate 

w + ek 



(9) 



(10) 



where Ak. 



Apgi/Jk and Ak,sc = Agcipk- Here we have 



added to Spg an additional piece Eq which is not accounted 
for by our (particle-particle) ladder diagrams. This leads to an 
"incoherent" background contribution which we will address 
later in the context of the cuprates. For the rest of the dis- 
cussion in the next two sections we set Eq = and, thereby, 
focus only on the superconducting and pseudogap terms. 
The spectral function can readily be computed from E = 



A(k, w) == -2 Im G(k, u + iO) , 



(11) 



which satisfies the sum rule f°° ^A('k,uj) — 1. We, thus, 
obtain a relatively simple expression for A(k, uj) which ap- 
plies below and above Tc, respectively: 
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FIG. 1. Behavior of the coherent contribution to the spectral function at various T, for wavevectors away from the nodes. 



A{k,u;) 
A(k,w) 



^^IpM^ + ^^f 



(12a) 



2A^7 



[uo^-ElY+^^iu-e^y 



(12b) 



From Eq. (12a), we see that the spectral function contains a 
zero at w — — ek below Tc, whereas it has no zero above Tc- 
This difference is responsible for the different thermodynam- 
ical behavior across Tc- 

In Fig. |I|, we plot the spectral function for ek = (on the 
Fermi surface) at different temperatures from slightly above 
Tc [Fig. Wa)] to temperatures within the superconducting 
phase [Fig. [l|(f)]. This figure is typical of situations in which 
there is a well established pseudogap. The figure can be 
viewed as representative of both s- and d-wave order param- 
eter symmetries. Hence the value of the wave-vector k is not 
particularly relevant, provided it is away from the nodal points 
in the d-wave case. For illustrative purposes, we take k at the 
anti-nodes, with ^{T) = Apg{Tc), and Apg{Tc) = 0.05 (in 
units of 4i||). In this way we ignore any T dependence in 7 
and, thus, single out the long range order effects associated 
with Asc- 

These figures give the first clear indications of the onset of 
"quasiparticle" coherence. Moreover, panel (a) helps to em- 
phasize an important component of our physical picture: the 
superconductor is not in a Fermi liquid state just above Tc, 
as can be seen by the non-Fermi liquid form for the spec- 
tral function. Just below Tc, the dramatic dip at 0.99Tc is a 
consequence of Bose condensation of q = pairs. Here, a 
very small condensate contribution nevertheless leads to the 
depletion of the spectral weight at the Fermi level, as shown 
in Fig. |l](b). As the temperature continues to decrease, and 
the superconducting gap increases, the two peaks in the spec- 
tral function become increasingly well separated, as plotted in 
Figs. n^(c)-(f). Even at the relatively high temperatures corre- 
sponding to T/Tc ~ 0.7, the spectral peaks are quite sharp 
— only slightly broadened relative to their BCS counterparts 
(where the spectral function is composed of two 6 functions). 



It should be stressed that lifetime effects via 7 do not lead 
to significant peak broadening. This follows from the fact that 
the imaginary part of the pseudogap self-energy at the peak 
location _Ek is given by (for ek = 0) 



7=7 



Av 



(^k 



kk|)^ + 7' 



A^ + 7^ 



(13) 



Since Eq. (ra) indicates that Apg vanishes as T ^ 0, the ef- 
fective peak width, determined by 7', decreases with decreas- 
ing T. It can be seen that below Tc, the spectral function in 



Eq. (12b), is very different from that obtained using a simple 
broadened BCS form; there is no true gap for the latter, in 
contrast to the present case. 

These spectral functions can be used to derive the density 
of states (per spin) as 



iVH = ^A(k,c.) 



(14) 



Moreover, it is expected that peak sharpening effects dis- 
cussed above for the spectral function are also reflected in the 
density of states. For siinplicity, we first consider the case of 
s-wave pairing. In Fig. H(a)-(f), the density of states is plot- 
ted for a quasi-2D s-wave superconductor, where the various 
energy gaps are taken to be the same as in Fig. [l|. Because 
of contributions from states with ek ^ 0, the narrow dips in 
Fig. [l|(b)-(c) do not show up here. However, as is evident, the 
density of states within the gap region decreases quickly, as 
the superconducting condensate develops. 

The rapid decrease of the density of states with decreasing 
T, in the vicinity of Tc, will be reflected in the behavior of the 
specific heat, Cy and, thereby, lead to the thermodynamic al 
signature of the phase transition. C„ may be obtained from 
Cy = dE/dT, where the energy E is calculated viaS 

E = 2TY, ^(*^n + 4 + M)G(k,ic^„) 

k,n 

= E/ ^ioJ + ek + 2^i)A{k,Lo)f{uj), (15) 

.. "'-00 ^'^ 
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FIG. 2. Effects of superconducting long range order on the behavior of the density of states as a function of temperature in a pseudogapped 
s-wave superconductor with n — 0.5. Here we take the parameters to be the same as those used in Fig. [l| At T/Tc ~ 0.7, as shown in (f), the 
density of states is close to that of strict BCS theory. 



where e^ = e^+fj. is the dispersion measured from the bottom 
of the band. It follows that 



0.3 



E ■ 



2^ 



[iu + ^,)N{u;) + K{L,)]f{L,), 



(16) 



where we have defined K{u!) = J2k ^k^C*, ^), which can be 
regarded as the contribution associated with the kinetic energy 
of the system. In this way, we obtain 



C„ 



dw 

2^ 



9/1 



NiL,)f{Lu) 



dT 

{uj + ^l)N{uJ) + K{uo) 
dN{oj) 



(uj + fi)- 



T 

dK{uj) 



dT 



dT 






(17) 



The first two terms on the right hand side lead to a "nor- 
mal metal-like" contribution to Cv/T which is proportional 
to N{uj) at low T. However, the third term arises because 
N[u)) depends on T. In this case, Cy/T no longer reflects 
the density of states. It is this term that will give rise to the 
specific heat discontinuity at T^. 

In Fig. we plot the temperature dependence of C^ in both 
(a) the weak coupling BCS case and (b) the moderate coupling 
pseudogap case with s-wave pairing. We choose, for definite- 
ness, the broadening ^{T) = T for the second of these calcu- 
lations. We also indicate in the insets, the respective tpmnpi-a- 
ture dependent excitation gaps, which have been assumedE3l£l 
in producing the figure. 

In both cases shown in Fig. g, the-specific heat jump arises 
from a discontinuity in dA^(aj)/dT,a associated with the on- 
set of superconducting order. However, for the BCS case, this 
derivative can be associated with a discontinuity in the deriva- 
tive of the excitation gap, via 



AC, 



BCS 



-N{Q) 



dA2 
dT" 



(18) 
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FIG. 3. Comparison of the temperature dependence of the spe- 
cific heat in the (a) weak coupling BCS case and (b) moderate cou- 
pling pseudogap case. Shown here are quasi-2D s-wave results, at 
n = 0.5, — gr/4f II = 0.5 and 0.6, respectively. The T dependence of 
the gap is shown as insets. 



By contrast, in the pseudogap case, the gap A and its deriva- 
tive dA/dT are presumed to be continuous across Tc as shown 
in the inset to Fig. ||(b), and in Fig. Qbelow, so that Eq. ( |l8| ) 
does not hold. Moreover, in this case, above but near T^ the 
temperature dependence in the density of states is still impor- 




T/T, 

FIG. 4. Temperature dependence of the excitation gaps for var- 
ious doping concentrations used for calculations in Fig. M. Here 
Ao.opt is the zero T gap at optimal doping x ~ 0.15. 



measurements for under- and optimally doped cupratesj2aEZl 
and places the van Hove singularity in a more correct position 
within the band. 

Finally, we turn to the phenomenological parameter 7 as 
well as to A. We presume that 7 changes from above to be- 
low Tc and in this way A (which is directly coupled to 7 via 
the set of coupled equations for G and T ) will have some, 
albeit small, structure in its temperature dependence at Tc, as 
seems to be the case experimentally. Our choice for the ex- 
citation gaps is shown in Fig. ^ and appears compatible with 
Figs. 8 and 9 in Ref. p2l As i^-oMKistent with scattering rate 
measurements in the literaprej2a E3 we take 7 oc T"^ below Tc 
and linear in T above TcE^ For the doping dependence, we 
assume that 7 varies inversely with A. This reflects the fact 
that when the gap is large, the available quasiparticle scatter- 
ing decreases. With these reasonable assumptions, along with 
the continuity of 7 at Tc, we obtain a simple form: 



7 



aT^/TcA , 
aTTc/A , 



(T < Tc) 
(T > Tc) 



(19) 



Here, the coefficient a < 1. 
adjustable parameter. 



This corresponds to our single 



tant due to the presence of an excitation gap above Tc- The lat-. 
ter leads to a decrease in diV(aj)/dT which is then reflectedEj 
in a decrease in Cy, slightly above Tc- At higher T, well away 
from Tc, where diV((jj)/dT tends gradually to zero, Cy is then 
controlled, as in a more typical "normal metal", by N{uj). 
We see, then, that the approach to the "normal metal" value 
is sharp for BCS, but because of the non-zero pseudogap, it 
is more gradual for case (b). An important consequence of 
these effects, is that the shape of the anomaly in Cy (shown in 
Fig. p](b)) is more characteristic of a A-like transition, although 
there is a precise step function discontinuity just at Tc- 



IV. APPLICATION TO THE CUPRATES 

The results obtained in Sec. [Il| are generally valid for both 
s- and d-wave cases, and can be readily applied to the d-wave 
cuprates. In this section we test the physical picture and the 
results obtained above, by studying the tunneling spectra and 
the specific heat behavior in the cuprates, as a function of dop- 
ing and of temperature. As in earlier workJijO we introduce a 
hole concentration dependence of the electronic energy scales 
by imposing the Mott constraint that the in-plane hopping in- 
tegral t|| (x) = tox, so that the plasma frequency vanishes as 
X ^ 0. As a result, the effective coupling strength —g/t\\{x) 
increases as the Mott insulator phase is approached. Here we 
assume (for simplicity) g{x) ^-.q and fit the one free parame- 
ter g/tQ to the phase diagram.ll3 

In order to compare with tunneling spectra, we introduce 
a slightly more realistic band structure which includes a next- 
nearest neighbor hopping term, — 2i'(l — cos k^ cos ky), in the 
band dispersion, ek, with i'/i|| « 0.4. This parameter choice 
gives rise to the hole-hke Fermi surface shape seen in ARPES 



A. Tunneling spectra 

Tunneling experiments were among the first to provide in- 
formation about the excitation gap — which measurements 
seem to be consistent with ARPES data.|j For a given density 
of states N{lo), the quasiparticle tunneling current across a 
superconductipc-insulator-normal (SIN) junction can be read- 
ily calculated,cd 



/; 



SIN 



= 2eNoT^ 



^N{uj)[f{cu-eV)-f{oj) 

2,71 



(20) 



where we have assumed a constant density of states, Nq, for 
the normal metal, and presumed that the tunneling matrix ele- 
ment To is isotropic. In reality, there may be some directional 
tunneling which will tend to accentuate the gap features, but 
we do not complicate our discussion here with these effects. 
At low T, one obtains 



d/ 
dF 



e'N.Ti 



N{eV) 



(21) 



SIN 



SO that the tunneling spectra and the density of states are 
equivalent, up to a multiplicative constant coefficient. At T 
comparable to Tc, however, the tunneling spectra reflect a 
more thermally broadened density of states. 

In Fig. g(a), we plot the SIN tunneling spectra, calculated 
for optimal doping {x « 0.15) at temperatures varying from 
above to below Tc. The van Hove singularity introduces a 
broad maximum in the spectra at high temperatures, as seen 
for the top curve in Fig. @(a). We see here thai (even for this 
optimal sample), as observed experimentally,^ the density of 
states contains (pseudo)gap like features which lead to two 
peaks. This is visible for temperatures well above Tc- A 
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FIG. 5. (a) Temperature and (b) doping dependence of tunneling spectra across an SIN junction. Shown in (a) are the dl/dV characteristics 
calculated for optimal doping at various temperatures from above to below T^. Shown in (b) are tunneling spectra at low T (at 0.2Tc) for 
various x. The units for dl/dV are e^A^oTo /4i||, where t\\ is evaluated at optimal doping. For clarity, the curves in (a) and (b) are vertically 
offset by 1.5 and 10, respectively. 



similar plot is presented in Fig. @(b), which shows at fixed 
low T = 0.2Tc how the spectrum evolves as a function of x. 
Both these plots appear in reasonable agreement wi|th| what is 
observed experimentally by Renner and coworkersc3 and by 
Miyakawa et a/B for Bi2Sr2CaCu208+^- (Bi2212). 



B. Specific lieat 

There is a substantial amcpnl of experimental data on the 
specific heat in the cuprates,oH although systematic studies 
come primarily from one experimental group.O We compare 
our numerical results with these data by plotting our calcula- 
tions for Cy/T in Fig. 0(a)-(f), from over- to underdoped sys- 
tems. As shown in these plots, the behavior of C^, is BCS-like 
in the overdoped regime. As the system passes from optimal 
doping towards underdoping, the behavior is more represen- 
tative of a A-like anomaly, as found in Fig. |](b) for the s-wave 
case. All these trends spem to be qualitatively consistent with 
experimental data.L3o In the underdoped regime at high T, 
we find a maximum in C^/T, near a temperature, T*, which 
may be associated with the onset of the pseudogap state [See, 
e.g. Fig. n(c)]. Finally, in contrast to Fig. H, at low T, the 
d-wave nodes lead to a larger quasiparticle specific heat, than 
for the s-wave case. 

The experimentally observed A-like anomaly of C^ at Tc 
has been interpreted previously as evidence for a Bose con- 
densation description.Q Here, in contrast, we see that within 
our generalized mean field theory, this anomaly naturally 



arises from the temperature dependence of the fermionic ex- 
citation gap which has some structure at, but persists above, 
Tc, as shown in Fig. 0. Thus, this is a property of supercon- 
ductors which have a well established pseudogap. That the 
experimental data, (which, except at extremely reduced hole 
concentrations), show a reasonably sharp (A-like) structure at 
Tc — seems to reinforce the general theme of this work — 
that corrections to BCS theory may be reasonably accounted 
for by an improved mean field theory, rather than by, say, in- 
cluding order parameter fluctuation effects. 



V. LOW T EXTRAPOLATION OF THE PSEUDOGAPPED 
NORMAL STATE 

A. Non-Fermi liquid based superconductivity 

The character of the extrapolated (T < Tc) "normal state" 
is at the core of many topical issues in high Tc superconduc- 
tivity. Understanding this state may shed light on the nature 
of Fermi liquid break-down above Tc- Moreover, the thermo- 
dynamics of this extrapolated phase provide a basis_£pr esti- 
mates of the condensation energy. This is deducedo by in- 
tegrating the difference between the entropy of the supercon- 
ducting state and that of the extrapolated normal state with 
respect to T. This extrapolated normal state also appears as 
a component of the free energy functional of conventional 
Landau-Ginzburg theory. Indeed, considerable attention has 
been paid recently to the condensation energy in the context 
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FIG. 6. Temperature dependence of the specific heat for various doping concentrations, calculated with a = 1/4 in Eq. h9l 



of determining the paid ng, "mechanism" in the high tempera- 
ture superconductors .EJCJ In this regard what is needed is the 
difference between the various "normal" and superconduct- 
ing state polarizabilities (e.g., magnetic and electric) which 
are thought to be responsible for the pairing. It is important 
to stress that the "normal state " used in computing the po- 
larizabilities should contain an excitation gap which is com=. 
patible with that found in the experimental data analysix^ 
with which the microscopically deduced condensation energy 
is compared. 

We have emphasized that within our physical picture there 
is an underlying (pseudo)gap in the normal state below Tc- 
A similar picture was independently deduced phenomenolog- 
ically from specific heat and magae.tic susceptibility measure- 
ments by Loram and co-workersE3 who arrived at equations 
like those of Section II [Eqs. (^ and (^]. However, they did 
not impose a self consistent condition on Apg as in Eq. (ra). 
Rather the quantity Apg (which they call Eg) is assumed to be 
T-independent. 

To expand on these issues we plot in Fig. ^ the calculated 
Cv/T and entropy 5* for (a)-(b) the BCS case, as compared 
with the counterparts obtained for the pseudogap supercon- 
ductor in (c) and (d). The dotted lines represent the Fermi 
liquid, i.e., linear extrapolation (FL). Figures pi a) and (b) re- 
affirm that this Fermi liquid extrapolation is sensible for the 
BCS case — Cy/T is a constant, and 5 is a straight line going 
through the origin. Panel (b) is useful in another respect: it 
shows how the entropy behaves as phase coherence is estab- 
lished. In general, the phase coherent state has a lower entropy 
than the extrapolated normal state. ■— ■ 

In the underdoped regime, Loram and co-workersEJ have 
stressed that entropy measurements lead one to infer that an 
excitation gap occurs above Tc- We analyze our calculated 
form of the entropy in a similar fashion. In contrast to the BCS 
case, for a pseudogap superconductor, the Fermi liquid extrap- 
olation of S is unphysical, approaching a negative value at 
low T, as shown by the dotted line in Fig. M(d). Here the solid 
and dotted lines separate around the temperature T* which we 



find to be around 1.5Tc. In order to obtain a thermodynami- 
cally consistent picture, then, the normal state must deviate 
from the FL line and this is accomplished by turning on an 
excitation gap at T < T*. 

The dashed lines in Fig. 0(c)-(d) show a more reasonable 
extrapolated normal state (labelled PG) which is equivalent to 
the solid line for T > Tc and distinct for T < Tc. This ex- 
trapolation is taken to be consistent with the conservation of 
entropy, S = L Cy/TdT, i.e., the shaded areas in (c) and 
(d). This construction for the T < Tc normal state, is similar 
to the procedure followed experimentally ,ej and, in effect, re- 
moves phase coherent contributions which enter via Age- This 
construction is by no means unique; all that is required is that 
the entropy of the extrapolated normal and of the supercon- 
ducting states be equal at T = Tc- As shown in the figure, we 
chose, for simplicity, a straight line extrapolation for Cy/T. 
Moreover, this choice is consistent with our expectation that 
there would be, in the "normal state", a finite intercept for 
Cy/T. 

Just as a gap is present above Tc, the underlying nor- 
mal phase below Tc (labelled PG) is to be distinguished 
from the FL extrapolation; it also contains an excitation 
gap. Indeed, this is consistent with what has been claimed 
experimentally :cJ a thermodynamically consistent picture of 
Cy reflects a gap in the T < Tc normal state spectrum, not di- 
rectly related to the condensate. While Figs. M(c)-(d) are sim- 
ilar to what is in the data in underdoped cuprates, in our anal- 
ysis the "normal" state Cy/T and 5* are linear and quadratic 
in T, respectively. Slightly different power laws have been as- 
sumed experimentally. A rough estimate of the condensation 
energy can be obtained from the integrated area between the 
solid line (for the superconducting state) and the dashed line 
(for the extrapolated normal state) in Fig. 0(d). It should also 
be noted that a more meaningful measure of the condensation 
energy is obtained by computing the magnetic field dependent 
Gibbs free energy. This is more complicated to implement 
both theoretically and experimentally. 
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B. Comparison with vortex core and Cv measurements 

The presence of a pseudogap in the underlying normal state 
of the superconducting phase is alsaconsistent with the obser- 
vations by Renner and co-workersS based on STM measure- 
ments within a vortex core. While one might be concerned 
about magnetic field, H, effects in interpreting these data, Jl 
should be noted that H appears to have a rather weak effectEll 
on pseudogap phenomena, as measured by T* and A(iJ). [In 
more overdoped samples the field dependence becomes more 
apparentCj]. By contrast Tc is more sensitive to H. 

Indeed, this weak dependence on H is often invoked in the 
literature as strong evidence against the "preformed" pair sce- 
nario. In the usual BCS case, pairs form precisely when phase 
coherence sets in. However, in the case of a pseudogap super- 
conductor where the coupling is stronger, pairs form above Tc 
without an underlying phase coherence. It is clear that a mag- 
netic field lowers Tc by destroying phase coherence. However, 
"preformed" pairs will survive H, leaving the excitation gap 
in tact. Stated alternatively, a magnetic field (just like mag- 
netic impurities) breaks time reversal symmetry, and there- 
fore makes it energetically unfavorable to form Cooper pairs 
which are comprised of time reversed single particle states. 
In contrast, finite momentum pair excitations, which are re- 
sponsible for the pseudogap in our approach, are not formed 
in time reversed states, and as such, they are not as susceptible 
to external magnetic fields or to magnetic impurities. 

In Fig. pi we plot our results for the SIN tunneling 
characteristics, dl/dV, and the computed entropy and spe- 
cific heat in_a|jp|Seudogap superconductor, and compare with 
experiment.EjQ'CJ To obtain the extrapolated normal state 



(called PG) in dl/dV, we set the superconducting order pa- 
rameter Asc to zero, but maintain the total excitation gap to 
be same as in a phase coherent, superconducting state — with 
non-zero Age (called SC). This procedure presumes that when 
the condensate is absent, the pseudogap Ap^ must correspond 
to the full excitation gap. Thus it should reflect the pairs which 
would otherwise be condensed. The characteristic behavior 
of dl/dV measured in an SIN configuration is presented as 
a comparison between theory (left) and experiment (right) in 
the upper panels of Fig. || Here the experimental curves are 
taken from Ref. g[ measured for underdoped Bi2212 inside 
(PG) and outside (SC) a vortex core, respectively. The non- 
Fermi liquid nature of the extrapolated normal state can be 
clearly seen. In a similar fashion, we show the comparisons 
for the extrapolated entropy and specific heat between our the- 
ory and experiments of Ref. p3 for Yo.8Cao.2Ba2Cu307_A-. 
Here 7 = Cy/T. The agreement between the theoretically 
computed curves and experimentally deduced curves provides 
reasonable support for the present theoretical picture. 



C. Discussion 

In this section we re-visit some of the issues raised in this 
paper and in experiments on Cy, tunneling, vortex core and 
related spectroscopies. In contrast to what has been presented 
up until now, here we are more qualitative and, in some in- 
stances, more speculative. j— . 

Intrinsic tunneling experimentsEj Considerable attention 
has been directed towards intrinsic tunneling experiments (in 
stacked layers), not only because they yield different results 



10 
^ 5 



1 ' 1 ' 1 ' 


1 ' 1 
-— sc 


f V 


PG 


\ /I 


1 , 1 



-0.4 -0.2 





y 



0.2 



0.4 




-300 -200 -100 100 200 300 

1 I 

Sample ' 



0.2 



fe 



O 



: 0.1 



CO 







T/Tr 



— sc 

■ -— PG 

; Cy/j/x 


1 ' 

x= 0.125 ; 


1 









50 



' ' ' ' r ■'•-■--■- 1 1 1 ' ' ' 


r •■-•■' ' 1 1 ' i ■ 1 


y M 


6=0.439 




E =220K : 

% 


--^^~- — 


■ !/'y^^ 


- 


1 1 1-1... , . 


t . . 1 



100 150 

T(K) 



200 



250 



FIG. 8. Extrapolated normal state (PG) and superconducting state (SC) contributions to SIN tunneling and thermodynamics (left), as well 
as comparispji with experiments (right) on tunneling for Bi2212 from Renner et a/.u and on specific heat for Yo.8Cao.2Ba2Cu307_i from 
Loram et alEj The theoretical SIN curve is calculated for T — Tc/2, while the experimental curves are measured outside (dashed) and inside 
(solid) a vortex core. 



from STMc3 EJ and from point contact/break junctioncj ex- 
periments, but also because they sometimes reveal an unex- 
pected sharp feature (or second peak in dl/dV) presumably 
associated with superconductivity. This peak occurs in addi- 
tion to a broader excitation gap feature (which is more like 
that found in single junction experiments) and it vanishes for 
T > Tc- There is, as yet, no complete convergence between 
different intrinsic Josephson junctions (IJJ) tunneling jexperi- 
ments. On overdoped samples Suzuki and co-workersEil have 
found anomalously sharp and large amplitude (second) peaks 
whose presence correlates with long range order, while, by 
contrast, Latyshevo and co-workers find only a single max- 
imum in dl/dV below Tc, which rather smoothly evolves 
into the normal state peak. However, for underdoped sam- 
ples, Krasnov and co-workersE3 report two maxima below Tc, 
with a much less pronounced sharp feature than in Ref. bTI 
(although the latter appears to be descended from the more 
anorrtaJous features reported earlier by Suzuki et al). Earlier 
worko by some of his same co-authors found only a single 
gap feature, *^was consistent with single junction data from 
other groups. Cju 

There is in all these experiments the possibility that the so- 
called superconducting peak is an artifact of self-heating or 
other non-equilibrium effects, which would be present when 
there is a non-zero critical current. Its absence in single junc- 
tion experiments would, otherwise, be difficult to explain. 
These latter experiments coiTelate well with ARPES. More- 
over, they also correlate with inferences from Cy and other 



bulk dataoEl We have no simple explanation for the sharp 
feature in tunneling. Within our approach there is a single 
excitation gap,_A above as well as below Tc- Although the 
critical currenttS Ic reflects the order parameter Age, which 
vanishes at Tc, this order parameter contribution is not ex- 
pected to show up in quasiparticle tunneling as a second gap 
feature. Of these IJJ experiments, the data which seems not 
incompatible with our pjicture is that of Yurgens et a/.£j and 
possibly Latyshev et al. ;l3 the latter authors, nevertheless, find 
much sharper maxima in dl/dV than we wouJd, have. 

Quantum, ^itical points: Both theoristsEa and experi- 
mentalistsEjO have recently turned their attention to quantum 
critical phase transitions. Moreover, these phenomena are as- 
sumed to bea-elated to pseudogap effects. Indeed, Loram and 
co-workersEil presume that Apg (which, in their approach, is 
taken to be temperature independent) is proportional to T*. 
Then at some critical doping concentration (x w 0.19), Apg 
appears to vanish and they infer that T* -^ 0. By contrast, we 
find Apg is more closely associated with (T* — Tc), which 
does not lead to a zero temperature phase transition, even 
when Apg vanishes. 

At low X, on the other hand, there may be something more 
dramatic like a first order or quantum critical phase transition 
going on — at the superconductor-insulator boundary. Here 
the excitation gap is maximum on one side of the bound- 
ary and yet the superconducting order parameter disappears 
on the other. [In the present picture this disappearaoce was 
found to arise from the localization of d-wave pairstj]. At 
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low X, this superconductor-insulatpi. transition also appears 
in the presence of a magnetic fielcO when the field is large 
enough to drive the system into the-normal phase. It also ap- 
pears with impurity pair breaking.EZi All three of these exper- 
iments may be interpreted as suggesting that the fermionic 
excitation gap A survives in the presence of pair breaking 
(by large fields or impurities), or low hole concentrations — 
thereby leading to an insulating fermionic excitation spec- 
trum. Some confirmation of this conjecture comes from NMR 
experimentscZl which seem to imply that A does not vary with 
H, once the pseudogap is well established. And the STM 
measurements,! in general, as well as inside a vortex core 
seem consistent with the observation that A is only weakly H 
dependent. One cannot, of course, ignore Mott insulating ef- 
fects at this superconductor-insulator boundary as well. How- 
ever, whatever physical mechanism is dominant, the fact that 
the superconductor-insulator transition is a robust feature as- 
sociated with the disappearance of superconductivity, whether 
by pair-breaking or by doping, needs to be addressed. 

Incoherent contributions to the spectral function in photoe- 
mission: . .. . 

Recent photoemission experimentsL3E3 indicate that the 
relative weight of the coherent contribution to the spectral 
function decreases rapidly with decreasing x. Ignoring the in- 
coherent term So in the self energy (as we have here) will nec- 
essarily affect any quantitative inferences about the systematic 
X dependence of the spectral function and, thereby, of Cy or 
the related condensation energy. This term can only be put in 
by hand in the present approach; it arises here from diagrams 
other than the particle-particle terms which give rise to the su- 
perconductivity and pseudogap. Moreover, the measured sys- 
tematic X dependence of the coherenL-snectral weight, which 
has been inferred from photoemissionj3E] is likely tpie con- 
sistent with the inferences based on thermodynamicsEi for the 
X dependence of the entropy 5* and related condensation en- 
ergy. However, when the contribution from Eq is sizeable, 
relative to the coherent terms, one cannot include it (by hand) 
without self-consistently also re-solving for the chemical po- 
tential /i and, hence, also for A and Tc, etc. This extensive 
numerical program would take us too far afield to implement 
here. 

In addition, this Eo contribution is needed to arrive at the 
well knowrO dip-hump features of photoemission. When 
this term is artificially added, we are able to obtain this lat- 
ter structure which will scale with the excitation gap A. In- 
deed, a dependence on A is plausible since we presume that 
So is associated with various (electron-hole) polarizabilities, 
which in the superconducting and pseudogap states reflect the 
non-vanishing excitation gap. It is, however, essential (to ob- 
tain dip-hump features) that the imaginary part of Eq turns on 
rather abruptly at frequencies somewhere between A and 2 A. 
Indeed, the step function model introduced in Ref. g^ seems 
to accomplish this quite well, but we know of no simple mi- 
croscopic mechanism which yields this rapid frequency onset. 



VI. CONCLUSIONS 

This paper has raised some issues which have a number 
of important ramifications. We suggest that specific heat and 
vortex core experiments have provided strong evidence that 
the normal state underlying the superconducting phase is not 
a Fermi liquid. Thus BCS theory cannot be applied to the 
underdoped cuprates, without some modification. Since we 
have, as yet, very little alternative to BCS, it is natural first 
to try to extend it slightly. We believe the simplest and most 
benign modification is to adopt the Leggett extension of the 
ground state wavefunction, Eq. (M, (which is applicable to 
weak and strong coupling g) and extend his scheme to finite 
T. Within this approach we werepabl£,in the past to perform a 
number of concrete calculationsJ13ilZl and, in the present pa- 
per, explore the behavior somewhat above and below Tc of the 
specific heat, Cy, and quasiparticle tunneling characteristics, 
dl/dV. This present study led us naturally to analyze the na- 
ture of superconducting phase coherence in the presence of a 
pseudogap. 

Our study of Cy and dl/dV is based on a Green's func- 
tion decoupling scheme chosen to be consistent with Eq. ([l|). 
The spectral functions which enter these two physical quan- 
tities depend, in turn, on the self energy which is the sum of 
Eqs. (m and dlQ), where throughout we have ignored the in- 
coherent term Eo, which enters through diagrams other than 
those responsible for the superconducting and pseudo-gaps. 
It should be noted that this form for E is different jfcpm that 
introduced phenomenologically by,Franz and MillisEa and by 
NormanEl One of these groups,E3 in particular, emphasized 
the effects of the temperature dependent scattering rate. Here, 
by contrast, we emphasize the effects of long range phase co- 
herence which sets in at Tc- 

Because of the breakdown of BCS theory, in a supercon- 
ductor with a pseudogap, the standard simplifications, such 
as Landau-ijinzburg expansions and Bogoliubov-de Gennes 
approachestj are not expected to hold, at least without some 
modifications. For the BCS case, the expansion in terms of a 
small order parameter which is identical to the excitation gap 
at Tc is possible. However, when the order parameter and the 
excitation gap are distinct as in the pseudogap case, there is 
no straightforward way to expand the free energy in terms of 
a small order parameter and to reflect the existence of a well 
established excitation gap simultaneously. 

It should be emphasized, finally, that, in the non-BCS su- 
perconductor, there is an important distinction between Tc and 
the zero temperature excitation gap A(0). In the strong cou- 
pling, but still feanipnic regime, as the pseudogap increases, 
Tc is suppressed.QElllj This observation allows us to respond 
to the widely repeatpi criticism of this "preformed" pair ap- 
proach: namely thatEll "the gap is not closely tied to the on- 
set of superconductivity" — as inferred from, say, the lack of 
magnetic field dependence in the former Here we claim that 
in contrast to BCS theory, the excitation gap A is expected 
to be robust with respect to standard pair-breaking perturba- 
tions (such as magnetic fields and impurity scattering), while 
the order parameter. Age, is not. As emphasized throughout 
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this paper, the distinction between these two quantities is an 
essential component of the "preformed" pair approach. 
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